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ABSTRACT 

We introduce a thermochemical kinetics and photochemical model. We use high-temperature bidi- 
rectional reaction rates for important H, C, O and N reactions (most importantly for CH4 to CO 
interconversion) , allowing us to attain thermochemical equilibrium, deep in an atmosphere, purely 
kinctically. This allows the chemical modeling of an entire atmosphere, from deep-atmosphere ther- 
mochemical equilibrium to the photochemically dominated regime. We use our model to explore the 
atmospheric chemistry of cooler (T e ff < 10 3 K) extrasolar giant planets. In particular, we choose to 
model the nearby hot Neptune GJ436b, the only planet in this temperature regime for which spectro- 
scopic measurements and estimates of chemical abundances now exist. Recent Spitzer measurements 
with retrieval have shown that methane is driven strongly out of equilibrium and is deeply depleted on 
the dayside of GJ 436b, whereas quenched carbon monoxide is abundant. This is surprising because 
GJ 436b is cooler than many of the heavily irradiated hot Jovians and thermally favorable for CH4, 
and thus requires an efficient mechanism for destroying it. We include realistic estimates of ultraviolet 
flux from the parent dM star GJ 436, to bound the direct photolysis and photosensitized depletion of 
CH 4 . While our models indicate fairly rich disequilibrium conditions are likely in cooler exoplanets 
over a range of planetary metallicities, we are unable to generate the conditions for substantial CH 4 
destruction. One possibility is an anomalous source of abundant H atoms between 0.01-1 bars (which 
attack CH4), but we cannot as yet identify an efficient means to produce these hot atoms. 

Subject headings: planetary systems — planets and satellites: atmospheres — planets and satellites: 
individual(GJ 436b) — methods: numerical — radiative transfer, disequilibrium 
chemistry 



1. INTRODUCTION 

Currently, transiting extrasolar planets offer virtu- 
ally exclusive [^] opportunities for observing physical and 
chemical states of exoplanetary atmospheres. Over the 
past four years, retrievals of atmospheric molecules from 
multicolor transit photometry (i.e. transit spectra) have 
compelled the development of progressively more sophis- 
ticated atmopheric models to interpret the observations 
and understand underlying chemical and dynamical pro- 
cesses. In particular, atmospheric-chemistry modeling is 
evolving from strictly thermo-equilibrium models with 

1 Correspondence to be directed to mrl@gps.caltech.edu or 
gv@s383.jpl.nasa.gov 

2 The exceptions to the exclusivity are the few young, self- 
luminous planets as in the HR8799 system 



stationary chemical species, to coupled models (Zahnle 
et al. 2009a,b; Line et al. 2010; Moses et al 2011) in- 
corporating thermo-kinetics, vertical transport, and pho- 
tochemistry. Thus far, such efforts have been devoted 
to hot- Jupiter planets, especially HD 209458b and HD 
189733b, due to their favorable transit depths and eclipse 
brightnesses and, therefore, far greater availability of 
observational data. However, with the recent retrieval 
of molecular abundances in the atmosphere of GJ 436b 
(Stevenson et al. 2010; Madhusudhan & Seager 2011), 
exoplanetary science is venturing into a new territory: 
hot-Neptune atmospheric chemistry. GJ 436b is bound 
to serve as a prototypical planet anchoring the theoret- 
ical framework for understanding the hot-Neptune class 
of exoplanets, much as how HD 209458b and HD 189733b 
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have for hot Jupiters. It is also the first planet with ob- 
servable thermal emission that transits an M star. M 
stars are of particular interest since they constitute the 
majority of stars in the solar neighborhood, and they 
have close-in habitable zones, which enhances radial- 
velocity detectability and transit observability; therefore, 
M stars present the best opportunities to discover and 
characterize rocky, potentially habitable exoplancts in 
the near future. GJ 436b and GJ 1214b provide the only 
present test cases for atmospheric chemistry of planets 
orbiting M dwarfs. Therefore, an era of intensive in- 
vestigations of this planet is commencing. This paper 
presents our application of a state-of-the-art model seam- 
lessly integrating thermo-kinetics, vertical transport, and 
photochemistry to simulate the atmospheric chemistry of 
GJ 436b in a similar manner to Visscher et al. (2010) 
and Moses et al. (2011) , along with realisitic estimates 
of UV fluxes for this planet. 

The first transiting hot Neptune discovered (Butler 
et al. 2004, Gillon et al. 2007), GJ 436b, revolves 
around an M dwarf merely 10 pc away from Earth and 
has received much attention due to its interesting or- 
bital dynamics (Ribas et al. 2008, Mardling 2008, Baty- 
gin et al. 2009), interior properties (Nettelmann et al. 
2010, Kramm et al. 2011), and atmospheric properties 
(Stevenson et al. 2010, Lewis et al. 2010, Madhusud- 
han & Seager 2011, Shabram et al. 2011). The slightly 
eccentric orbit (eccentricity = 0.16) has a mean orbital 
radius of 0.0287 AU (Torres et al. 2008), and the planet 
probably has a pseudo-synchronous rotation (Deming et 
al. 2007). The planet's mass is 23 M®, and its den- 
sity of 1.7 g/cm 3 resembles that of the ice-giant Neptune 
(1.63 g/cm 3 ). Analyses of its mass-radius relationship 
and transit depth indicates a layer of H/He dominated 
atmosphere is clearly required (Figucira et al. 2009; Net- 
telmann et al. 2010; Rogers & Seager 2010). The host 
star has an effective temperature of ~ 3400 K and an 
estimated age of 3 - 9 Gyr (Torres et al. 2008). Assum- 
ing zero albedo and global thermal re-distribution, the 
planet's effective temperature is 650 K. Of the confirmed 
transiting exoplanets (Wright et al. 2011), GJ 436b is 
one of the least irradiated and has one of the coolest 
atmospheres. Therefore, this planet represents a signifi- 
cant departure from hot Jupiters in terms of size, thermal 
environment, and UV flux. 

Although GJ 436b was discovered in 2004 (Butler, by 
radial velocity), it was not until 2010 that a retrieval 
of explicit molecular abundances in its atmosphere was 
reported (Stevenson et al. 2010), where six channels 
of secondary-eclipse photometry data ranging from 3.6 
to 24 /jm were analyzed by generating ~ 10 6 simulated 
spectra using varying combinations of molecular com- 
positions and temperature profiles to find the best fit 
to observations. A more recent paper (Madhusudhan & 
Seager 2011) provides further details and updated re- 
sults of a re-analysis of the same dataset using the same 
general retrieval method. In short, 10 6 combinations 
of ten physio-chemical free parameters, each spanning a 
large range of values, were used to generate synthetic 
dayside-emission spectra. In each of the 10 6 scenar- 
ios, six of the ten parameters were used to define the 
temperature-pressure (T-P) profile, whereas the other 
four parameters specified vertically uniform abundances 



of four molecules: H2O, CO, CH4, and CO2. Addi- 
tionally, the 1-D atmospheric model restricted the ra- 
tio of emergent flux output to incident stellar flux input 
on the day side to within the range between zero and 
unity. Given six data points and ten free paramters, the 
retrieval problem was mathematically undcrdctcrmincd. 
Nonetheless, sampling a million points in parameter- 
phase space allowed the authors to examine the joint 
probability contours, as defined by the goodness-of-fit 
(chi-square) function, projected on multiple-parameter 
spaces. Furthermore, by placing physical-plausibility 
constraints (in consideration of believable departures 
from thcrmo-cquilibrium chemistry) on the molecular 
abundances, the authors were able to confine the physical 
space to a fairly narrow, "best-fit," range for chi-square 
< 3. Depending on the wavelength, the photospheric al- 
titude varies from 9 bar to 0.2 bar levels. The main con- 
clusions are as follows: 1) temperature inversion is ruled 
out (i.e. no stratosphere); 2) 6 ppm (parts per million) 
is the absolute upper limit for CH4 abundance; 3) 300 
ppm is the absolute upper limit for H2O abundance; 4) 
CO2 and CO abundances are anti-correlated; 5) taking 
physical-plausibility into consideration, the best-fit spec- 
trum represents Ah 2 o = 100 ppm, Ach 4 = 1 ppni, Xqq 
= 7000 ppm, and Xqo 2 = 6 ppm, where A; is the number 
density of molecule i divided by that of H2. Also, note 
that even in the best-fit scenario, Xco 2 can range any- 
where from 1 - 100 ppm. The Stevenson et al. (2010) 
and the Madhusudhan & Seager (2011) efforts are the 
most comprehensive studies of atmospheric composition 
on GJ 436b thus far. 

From a theoretical point of view, the preceding abun- 
dance limits and values pose a very interesting challenge 
due to their drastic departures from thermo-equilibrium 
predictions, which indicate the following rough-order-of- 
magnitude values: Jh 2 o = 1000 (3xl0 4 ) ppm, Xqh 4 = 
1000 (10 4 ) ppm, A co = 60 (10 4 ) ppm, and X C o 2 = 0.1 
(1000) ppm for lx (50x) solar metallicities at <~ lbar. 
In either metallicity scenario, water and methane remain 
abundant (> 1000 ppm), whereas the retrieval shows wa- 
ter being relatively depleted and methane being drasti- 
cally depleted. Moreover, the thermo-equilibrium abun- 
dances of carbon monoxide and carbon dioxide are posi- 
tively correlated (either both low in the 1 x case or both 
high in the 50 x case), in contrast with the retrieval's 
anti-correlation. In particular, the retrieved partioning 
of carbon overwhelmingly in oxidized species amidst a 
hydrogen-dominated (reducing), temperate atmosphere 
is very surprising. For instance, at 1-bar pressure and 
solar metallicity, CH4 is the thermodynamically domi- 
nant carbon-bearing molecule for temperatures less than 
1100 K (Lodders & Fegley 2002). The common practices 
of simply adjusting metallicity and/or the C/O ratio can- 
not simultaneously reconcile these discrepancies. There- 
fore, one must investigate disequilibrium mechanisms. 

Madhusudhan & Seager (2011) posited that high 
metallicity combined with vertical mixing can explain 
the disequilibrium abundance of carbon oxides. Basi- 
cally, enhanced metallicity (<~ 10 x solar) can provide 
the requisite abundance of CO2. Since equilibrium CO 
abundance drops sharply with respect to temperature 
(Lodders & Fegley 2002) the retrieved uniformly high 
abundance of CO requires eddy mixing to populate up- 
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per, cooler, atmospheric layers. However, vertical eddy 
mixing alone cannot explain the large depletion of CH4 
due to its innately high thermochemical abundance in 
the deep atmosphere. Therefore, Madhusudhan & Seager 
(2011) invoked photochemistry as the potential culprit, 
based on Zahnle et al.'s (2009a, b) studies of photochem- 
istry on hot Jupiters. In such a scheme, photosensitized 
sulfur chemistry produces atomic H, which then destroys 
CH4 to form higher hydrocarbons. However, the Zahnle 
et al. (2009a, b) model uses solar-type stellar irradiancc 
and an isothermal atmsophere (i.e. constant tempera- 
ture versus altitude). As such, neither the photochemi- 
cal driver nor the thermal environment is tailored for our 
planet in question. More severely, Moses et al. (2011) 
pointed out that a typo in a key rate coefficient in the 
Zahnle et al. (2009a,b) model caused the apparent con- 
version of methane into higher hydrocarbons at pressures 
larger than 1 mbar. Generally speaking, at pressures 
larger than 1 mbar in a hydrogen-abundant atmosphere, 
hydrogenation of unsaturated hydrocarbons and reaction 
intermediates efficiently recycle species back to methane, 
preventing its large- scale destruction. Moses et al. 
(2011) also discussed the inadequacies of isothermal at- 
mospheric models due to their suppression of transport- 
induced quenching. Hence, the observed CH4 depletion 
still awaits adequate explanation. The low abundance of 
H 2 also has not been addressed. 

In addition to secondary eclipse observations, primary 
transit observations of GJ 436b exist as well (Pont et al. 
2009, Ballard et al. 2010, Beaulieu et al. 2011, Knutson 
et al 2011), and various groups have analyzed them to 
retrieve molecular abundances in the planet's termina- 
tor regions (Beaulieu et al 2011, Knutson et al 2011). In 
contrast to the secondary-eclipse retrieval, Beaulieu et al. 
(2011) were able to fit a compendium of their and Ballard 
et al.'s tansit observations between 0.5 and 9 ^.m with 
500 ppm CH4 in a H2 atmosphere, and finding no clear 
evidence for CO or CO2. Moreover, Beaulieu et al. pre- 
sented that a methane-rich atmosphere, with tempera- 
ture inversion, can be consistent with the said secondary- 
eclipse data as well (but see Shabram et al. 2011). More 
recently, Knutson et al. acquired Spitzer transit pho- 
tometry at 3.6, 4.5, and 8.0 /im during 11 visits. The 
multiple- visit data showed high transit-depth variability, 
which the authors attribute to potential stellar activity in 
the dM host. They did not find any compelling evidence 
for methane, and data excluding ones believed to be most 
affected by stellar activity appear to place an upper limit 
of 10 ppm for methane mixing ratio. The best-fit spec- 
trum to this select data set assumes 1000 ppm H2O, 1000 
ppm CO, 1 ppm CH4, with CO2 abundance poorly con- 
strained, roughly in agreement with Madhusudhan et al. 
Therefore, primary-transit data is currently inconclusive 
due to different interpretations by different groups. 

Our primary goal is to advance the fundamental un- 
derstanding of processes impacting the chemical state of 
GJ 436b by developing a 1-D atmospheric model that 
integrates all of the aforementioned equilibrium and dis- 
equilibrium processes. An important aspect of our model 
is the seamless integration of thermochemistry, kinetics, 
vertical mixing, and photochemistry in a manner that 
directly follows from Visscher et al. (2010), and con- 
temporaneously with Moses et al. (2011), obviating the 
conventional quench-level estimation (Prinn & Barshay 



1977). 

The quench-level approach assumes that the deep at- 
mosphere is in thermochemical equilibrium because high 
temperatures provide sufficient kinetic energy to over- 
come reaction barriers in either direction. However, 
as vertical tranport lifts a gas parcel to cooler, higher 
altitudes, chemistry becomes rate limited rather than 
thcrmodynamically determined. There comes a point 
in altitude where the kinetic conversion time scale be- 
comes slower than the transport time scale, and the rate- 
limiting reaction for a molecule of interest is not allowed 
time to reach completion. At altitudes above this point, 
the molecule's concentration is frozen/quenched (there- 
fore, the term "quench level"). In effect, the quench- 
level approach partitions the atmosphere into two parts: 
below the quench level, thermochemical equilibrium de- 
termines chemical abundances; above the quench level, 
molecular abundances are uniform versus altitude, with 
values equal to the equilibrium value at the appropriate 
quench level for each species. Although this approach 
has a long record of success (e.g., Prinn & Barshay 1977; 
Smith 1998; Griffith & Yelle 1999; Saumon et al. 2003; 
2006; 2007; Hubeny & Burrows 2007; Cooper & Show- 
man 2006), it does have some limiting assumptions and 
caveats that require great judiciousness. Specifically, one 
needs to determine the appropriate rate-limiting, inter- 
conversion reaction for each set of coupled species of in- 
terest (e.g., interconversion between CH4 & CO). The 
correct reaction choice is not always readily apparent (see 
e.g., Visscher et al. 2010) and the appropriate length 
scale for deriving the mixing time scale from the verti- 
cal eddy diffusion coefficient (K zz ) is still under some 
debate. Furthermore, since a basic assumption is that 
temperature decreases with altitude, atmospheric tem- 
perature inversions can complicate matters. 

Therefore, we implemented a fully reversible kinetic 
model in the following manner. Every measured forward 
reaction rate in our list is reversed using the equilibrium 
constant and the principle of microscopic reversibility. 
Given enough pathways, both forward and backwards, 
a given set of chemical species will reach thermochem- 
ical equilibrium, kinctically. This provides a seamless 
transition from the thermochemical equilibrium regime 
to the disequilibrium-dominated regimes. We can inves- 
tigate the disequilibrium effects on atmospheric compo- 
sition in a much more holistic, systematic manner, com- 
pared to heuristically identifying plausible disequilibrium 
processes. 

In the remainder of this manuscript we describe the dis- 
equilibrium processes that may be occurring in GJ436b's 
atmosphere. In §2 we describe thermochemical and 
chemical kinetics models as well as our estimate for the 
stellar UV flux. In §3 we show the modeling results as 
well as a description of the important reaction schemes 
governing the abundances of various species. Finally in 
§4 we discuss the relevant implications and conclude. 

2. DESCRIPTION OF MODELS 

We use joint thermochemistry and "1-D chemical- 
kinetics with photochemistry" models to study the at- 
mosphere's departure from thermal equilibrium. Exter- 
nal inputs to our models are the metals fraction (de- 
noted further on by (), the pressure and temperature 
(T-P) profile, the eddy diffusion coefficient profile, and 
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the incident stellar flux; note that we fix the T-P pro- 
file and the chemistry is decoupled from it, i.e., there 
is no self-consistent, radiative-convective adjustment of 
temperature structure when the chemistry is evolved to- 
wards steady state. We initialize the 1-D atmospheres 
using the NASA Chemical Equilibrium with Applications 
(CEA) model (Gordon & McBride 1996). Given the ini- 
tial elemental abundances of H, He, C, O, N, and S in an 
atmospheric layer, along with the layer's pressure and 
temperature, CEA uses a Gibbs free-energy minimiza- 
tion and mass balance routine to calculate the equilib- 
rium species abundances. 

Whereas chemical equilibrium concentrations are use- 
ful for initializing the atmosphere, they do not pro- 
vide the correct chemical state above pressure levels of 
~ 10 bars (Prinn & Barshay 1977; Griffith & Yelle 1999; 
Cooper & Showman 2006; Line et al. 2010; Moses et al. 
2011). We simply supply the equilibirum mixing-ratios 
as boundary conditions in the deep atmosphere for the 
kinetics calculations, and thereafter evolve the chemi- 
cal state over multiple timesteps until a steady state is 
reached. 

The computations are carried out with the Cal- 
tech/JPL photochemical and kinetics model, KINETICS 
(a fully implicit, finite difference code), which solves the 
coupled continuity equations for each involved species, 
and includes transport via molecular and eddy diffusion 
(Allen et al. 1981; Yung et al. 1984; Gladstone et al. 
1996; Moses et al. 2005). We use the H, C, and O chemi- 
cal reaction list originally described in Liang et al.(2003; 
2004) and references therein updated to high tempera- 
tures, recently augmented with a set of N reactions. We 
have not included the chemistry of sulfur in any great 
detail, because much of its kinetics is poorly constrained 
(see e.g., Moses et al. 1996). However we do consider 
a small, but well measured, set of H 2 S reactions. This 
helps us appraise if and how the introduction of S affects 
the abundances of the main molecular reservoirs of H, C, 
N, O such as CH 4 . 

We use high temperature rate coefficients for reac- 
tions from Line et al. (2010). All reactions are bidi- 
rectional, and we reverse them by calculating the back- 
reaction rates using thermodynamic data (see Table SI). 
With appropriate reaction pathways and proper rates for 
the back-reactions, the models can converge to chemical 
equilibrium purely kinctically in the deep planetary at- 
mosphere where reaction timescales are short compared 
to transport timescales, and photochemical reactions are 
unimportant. As mentioned earlier, this removes the 
cumbersome requirement of having to choose a lower 
boundary for individual species through ad hoc quench- 
level arguments (Prinn & Barshay 1977; Smith et al. 
1998). 

We solve for 51 hydrogen, carbon, oxygen and nitrogen 
bearing species including H, He, H 2 , C, CH, 1 CH 2 , 3 CH 2 , 
CH 3 , CH 4 , C 2 , C 2 H, C 2 H 2 , C 2 H 3 , C 2 H4, C 2 H 5 , C 2 H 6 , O, 
O(ID), 2 , OH, H 2 0, CO, C02, HCO, H 2 CO, CH 2 OH, 
CH 3 0, CH3OH, HCCO, H 2 CCO, CH3CO, CH3CHO, 
C 2 H 4 OH, N, N 2 , NH, NH 2 , NH 3 , N 2 H, N 2 H 2 , N 2 H 3 , 
N 2 H 4 , NO, HNO, NCO, HCN, CN, CH 3 NH 2 , CH 2 NH 2 , 
CH 2 NH, H 2 CN, with a total of - 700 reactions, 55 of 
which are photolysis reactions. The chemical pathway 
for reducing CO to CH 4 , described recently for Jupiter's 
deep atmosphere (Visscher et al. 2010), is included in our 
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Fig. 1. — Estimated temperature profiles for GJ 436b. The 
dashed profile is the disk averaged dayside profile retrieved by 
Stevenson et al. (2010). The solid curve is the f = 1 profile from 
the global circulation model of Lewis et al. (2010). We use the 
latter T-P profile for our chemical models. 

reaction list, along with the reverse pathways for CH 4 to 
CO oxidation. Photolysis absorption cross sections are 
from Moses et al. (2005) and the thermodynamic data 
(i.e. the compilation of entropies and enthalpies) used to 
reverse the kinetic rate coefficients are from JANAF and 
CEA thermobuild databases; e.g., CEA uses data from 
Chase et al. (1998) and Gurvich et al. (1989) (see Zehe 
et al. 2002). 

2.1. Model Parameters 

We model a large pressure and altitude range, 10 3 to 
lO" 11 bars (-5000 km or -0.2 R p from the 1 bar level), 
so as to capture the three major atmospheric regimes 
and the transitions between them. These three domi- 
nant portions of the atmosphere are - the thermal equi- 
librium regime in the deep hot atmosphere, the eddy 
transport dominated regime at intermediate pressures, 
and the photochemical regime at low pressures. A total 
of 190 pressure levels, uniform in logarithmic space, are 
used between the abovementioned levels, giving a reso- 
lution of about 14 levels per decade of pressure. Alti- 
tudes above the homopause remain relatively cool in our 
models, and we disregard the possibility of a hot thermo- 
sphere despite the models extend up to exophere levels 
at 10 -11 bars; this simplification has little or no bear- 
ing on the state of the atmosphere below the homopause 
(P — 1/ibar). We adopt the ( — 1 T-P profile from Lewis 
et al. (2010) (see Figure 1), noting its similarity to the 
T-P profile retrieved in Madhusudhan & Seager (2011) 
and Stevenson et al. (2010). Whereas GJ 436 itself is 
slightly subsolar in abundances (Bean et al. 2006), we 
allow for a span of planetary metallicities, covering the 
cases C = 0.1, 1, 50, and allowing for the possibility that 
the planet is either enriched or depleted; we used solar 
abundances from the standard text of Yung & DeMore 
(1999fl For non-solar atmospheres we tune the frac- 

2 Yung & DeMore (1999) tabulate the abundances of Anders & 
Ebihara (1982). These values predate the more recent downward 
revision of elements C, O etc. in the Solar photosphere (reviewed 
in Asplund et al. 2009). Our C/H, O/H, N/H and S/H ratios are 
a factor 1.66, 1.52, 1.35 and 1.43 higher than those recommended 
in Asplund et al. (2009). On this revised scale we are modeling a 
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tions of C, N, O, and S relative to H but not relative to 
each-other (e.g., C/O, N/O, S/O, are always fixed). 

The eddy diffusion strength (parameterized by a co- 
effcient, K zz ) determines the pressure level at which a 
species is chemically quenched. At the quench level for 
chemical X, the timescale for vertical transport (T trans ) 
equals the chemical loss timescale (T c hem,x)- Above that 
level, which includes the visible portion of the atmo- 
sphere, the mixing virtually "freezes" the concentration 
of that species. Below the quench level, T c hem,x "C 
T trans-, an d thermochemical balance is achieved. Line et 
al. (2010) and Moses et al. (2011) have used piecewise 
estimates of the eddy diffusion profiles, K ZZ (P). The 
recipe has been to estimate K zz in the deep adiabatic tro- 
posphere (~ 10 3 bars) using mixing length theories (e.g., 
Flasar & Gierasch 1977) and stitch this to global circu- 
lation model (GCM) derived profiles obtained by multi- 
plying the (horizontally averaged) GCM vertical winds of 
Showman et al. (2009) by the local scale height. Lewis et 
al. (2010) apply this procedure to their GJ 436b circula- 
tion model, and estimate that K zz increases from ~ 10 8 
at depth (100 bars) to 10 11 cm 2 s~ x at lower pressures 
(1 mbar). 

Such procedures have gnawing uncertainties - for ex- 
ample, the appropriate eddy mixing length may only be a 
fraction of the scale height, or the vertical wind strengths 
could well be overestimated. Smith (1998) has demon- 
strated theoretically that using an eddy length scale 
equal to the scale height is inappropriate, and may lead 
to gross over-estimates of the length scale (L) and the 
timescale (rt rans = L 2 /K zz ). Herein, we simplify mat- 
ters by choosing a constant K ZZ (P) = 10 8 cm 2 s _1 pro- 
file; this value is similar to that for the deep atmosphere 
in the Lewis et al. GCM. This simplification has a couple 
of redeeming features. First, this gives quench levels sim- 
ilar to those that would be derived had we used a GCM- 
inspired K zz profile. Second, whereas a low K zz may 
underestimate the mixing strength at higher altitudes, it 
has the effect more lethargic replenishment of methane 
and other photodissociated species from the lower atmo- 
sphere (it bolsters the photochemical timescale, relative 

tO T£ rans ). 

2.2. The Ultraviolet Emission from GJ 436 

dM stars such as GJ 436 show very little photospheric 
emission in the near to far ultraviolet (UV). Neverthe- 
less, non-radiative energetic processes can transport en- 
ergy to power a hot outer atmosphere, and this energy is 
partially dissipated in the form of cooling, chromospheric 
UV emission. Because the UV emission levels depend on 
many factors, ab initio estimates of it are difficult. We 
use GALEX and ROSAT derived estimates for GJ 436 
and combine these with a T e f f ~ 3400 K continuum from 
the stellar photosphere. This combined emission is used 
to drive photochemical reactions in GJ 436b. 

In the planetary atmosphere both H2 and He are 
weak absorbers relative to other molecular species, but 
are enormously more abundant. Helium ceases to ab- 
sorb longwards of 500 A, and H2 longwards of 1000 A. 
Methane, a carbon reservoir and the molecule of partic- 
ular interest herein, has a large absorption cross-section 

planet with £ ~ 0.16, 1.6, 80. This was brought to our attention by 
the anonymous referee. 



shortwards of 1600 A. Whereas methane (and water) is 
largely shielded by H2 and He from very shortwave radi- 
ation, it is photodissociated by radiation between 1000- 
1600 A, and is therefore susceptible to possible intense H 
I Ly a (A = 1216 A) from the M star host. Longwards of 
A = 1600 A, direct photolysis of methane dwindles due to 
a combination of the falling cross-section and weak stellar 
flux. Hydrogen sulfide photodissociates at much longer 
wavelengths, A < 2600 A, and if present in substantial 
quantities, is poorly shielded by other reservoir molecules 
H 2 , CH4, H 2 0, etc. H 2 S photolysis and the resultant hot 
atomic hydrogen may be influential if A ~ 2600 A pho- 
tons can penetrate deep into the planetary atmosphere 
(more in §3.3.5). 

GJ 436 is detected in a GALEX survey exposure in 
the near UV channel with flux f nuv = 21.0 ± 3.7 //Jy 
(near-UV channel, A = 2267 A, A\ FWHM = 616 A). It 
is undetected in the GALEX far UV band, with a 3er 
upper limit of ff uv < 24 /iJy (far-UV channel, A = 1516 
A, AX FWHM = 270 A). These can be converted to inci- 
dent UV photon fluxes at the mean orbital separation of 
GJ 436b. The near UV detection implies a flux of 9 x 10 10 
photons cur 2 s" 1 A" 1 A = 1960 - 2580 A at the plan- 
etary substellar point. This dosage at GJ 436b is about 
0.2 PELs (present-Earth-levels) ; mean Solar photon flux 
at Earth is 4.7 x 10 11 photons cm~ 2 s _1 A^ 1 between 
2000-2500 A (Yung & DeMore 1999). The 3cr flux upper 
bound (GALEX far-UV channel) is < 1.3 x 10 11 photons 
cm~ 2 s _1 A -1 A = 1450 - 1650A; this is just a factor 
of two higher than present-Earth-levels in an equivalent 
passband. 

H Ly a emission can be powerful in the upper chro- 
mospheres of cool stars. Because it is strongly ab- 
sorbed in the interstellar medium, direct line strength 
estimates are difficult. We make an indirect determi- 
nation based on empirical correlations with soft X-ray 
fluxes. Soft X-ray emission from GJ 436 has been ob- 
served in the Rosat All Sky Survey (Hiinsch et al. 1999), 
with f x ~ 5.4 x 10" 14 erg cur 2 s" 1 (0.1-2.4 keV; 
Rosat PSPC), implying a fractional X-ray luminosity of 
Lx/Lboi ~ 8 x 10 -6 ; this fraction is a factor ~ 100 lower 
than that observed from the most active dM stars and is 
consistent with GJ436b's estimated advanced age, 6 ± 3 
Gyr. More recent XMM-Newton EPIC measurements 
(Sanz-Forcada et al. 2010) give a factor of 8 lower L x , 
which may well be due to X-ray activity. Herein, we 
adopt the ROSAT flux because larger X-ray fluxes imply 
proportionally larger Ly a fluxes. 

To estimate the Ly a output, we use an an empirical 
correlation of the X-ray and Ly a emission of stars, de- 
rived from stellar samples that include several late type 
stars (e.g., Landsman & Simon 1993 and Woods et al. 
2004; in these papers, measurements of Ly a lines were 
made from International Ultraviolet Explorer and Hub- 
ble Space Telescope spectra, after applying a model based 
correction of ISM absorption) . Inverting the Woods et al. 
(2004) empirical power law, \ogF x ~ 2.2 log Fi ya — 7.76, 
we determine a photon flux of fhya ~ 1.5 x 10 14 pho- 
tons cm" 2 s -1 at GJ 43610 The Solar H Ly a flux at 

3 Very recently, Ehrenreich et al. (2011) estimate a Ly a flux 
using HST-STIS observations of GJ436. Their estimated line flux 
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Earth is ~ 10 12 photons cm -2 s _1 , a factor 100 lower. 
The reliability of X-ray derived Lyman a line flux may 
be assessed by comparing F^ ya with the GJ 436b's H 
a line flux. H a observed in GJ436 in absorption, with 
an equivalent width of 0.32 A (Palomar-Michigan State 
Nearby Star Spectroscopic Survey; Gizis, Reid & Haw- 
ley 2002), implies a line flux of Fn a — 2 x 10 5 erg cm -2 
s _1 , and a line strength ratio of H Ly a to H a of 2.2. 
For dM stars, where H Ly a is seen in emission and for 
which the intrinsic Ly a line strengths have been mea- 
sured, this line strength ratio varies between 3-5, with 
some stars having ratios as low as 2 and others as high 
as 8 (Doyle et al. 1997). 

3. CHEMICAL MODEL RESULTS 

3.1. Thermochemical Equilibrium 

Equilibrium vertical mixing ratios for the three metal- 
licity cases are shown in Figure 2: these are sub-solar 
( = 0.1, solar C = 1 an d super-solar £ = 50 heavy elemen- 
tal abundances. Because GJ 436b is significantly cooler 
than HD 189733b and HD 209458b, CH 4 is the thermo- 
chemically favored carbon carrier; higher effective tem- 
peratures drive equilibrium towards CO in the two hot 
Jupitcrs. The thermochemical abundances of CH4, CO 
and H2O along the T-P profile are readily understood 
through the net reaction 



CO + 3H 2 ^ CH 4 + H 2 
along with the Law of Mass Action: 
Xch4,Xh 2 o 1 



X CO X 3 H2 P 2 



K eq {T) 



(1) 



(2) 



derived by minimizing the Gibbs free energy of net reac- 
tion in (1), with the mixing ratio Xi of species i, with 
ambient pressure P, and a temperature dependent equi- 
librium constant K eq (T); the T dependence is governed 
by the van 't Hoff equation (AG = — RT log K eq , with 
AG as the standard Gibbs free energy change). At a 
given pressure P, K eq (T) behaves in a manner that rising 
T drives the equilibrium towards CO. At a fixed T, in- 
creasing/decreasing pressures favor higher CH4/CO con- 
centrations. These relationships are exemplified in the 
C, = 1 equilibrium profiles shown in Figure 2 (middle 
panel). As P and T decrease along the adiabat between 
1000— 100 bars, the equilibrium constant dominates over 
the adverse P 2 dependence, resulting in a drop in the CO 
fraction. In the isothermal region between 10—1 bars, 
decreasing pressure now favors the production of CO. Be- 
tween 1 bar and « 10~ 2 bars, the CO fraction falls be- 
cause of the rapid decrease in temperature with altitude. 
At levels above the ~ 10~ 2 level the temperature struc- 
ture is nearly isothermal, and the decreasing pressure 
favors higher CO fractions. Similarly, NH3 is the favored 
N carrier deep in the atmosphere, but is less favored at 
lower atmospheric pressures. Sulfur can be predominant 
as H 2 S, HS, or S depending on pressure and temperature, 
but for conditions prevalent in GJ 436b, gas phase H 2 S 
is the dominant sulfur reservoir and its concentration is 
unaffected by the temperature structure. Heavier hy- 
drocarbons, such as ethane (C 2 Hg), are relatively scarce 

is a factor 1.5 X smaller than the estimate based on h x used herein 
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Fig. 2. — Thermochemical equilibrium vertical distributions for 
abundant H, C, O, N, and S species assuming the temperature pro- 
file in Figure 2. Three metallicity cases are shown (f = 0.1, 1, and 
50, from top to bottom). The thermochemical equilibrium mixing 
ratios are derived using the CEA Gibbs free energy minimization 
code for each atmospheric T-P level. 



any pressure or temperature (but more common at the 
highest metallicities) . 

Enriching the atmosphere to £ = 50 increases the mix- 
ing ratios of the reservoir species in proportion, however 
the shapes of the vertical profiles are much the same as 
for solar metallicities. Similarly, decreasing the metallic- 
ity of the atmosphere to £ = 0.1 lowers the mixing ratios 
of the heavy gases, by a factor ~ £ for CH4 and £ 2 for CO 
etc. The shapes of vertical distributions are nonetheless 
preserved, and relatively insensitive to £. 

For all three metallicity cases considered, the chemi- 
cal equilibrium abundances of CH4 and H 2 stay rela- 
tively high there is always enough hydrogen present 
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to build these molecules. One can imagine an extreme 
situation where H is highly depleted, but such an atmo- 
sphere would be incompatible with the observed plane- 
tary radius. Conversely, the planet could be impover- 
ished in metals to greatly subsolar levels ( -C 0.1, al- 
though unreasonably low metallicities (< 1 x 10~ 5 x so- 
lar) would be required to deplete CH4 and other common 
molecules to levels below 1 ppm. These simple cases serve 
to show that, based solely on chemical thermodynamics, 
CH 4 has to be relatively abundant in GJ 436b and other 
T eff = 500 - 1000 K H-rich planets. 

3.2. Vertical Mixing & Chemical Quenching 

Vertical turbulent mixing has been invoked to explain 
the anomalously large observed abundance of CO in 
Jupiter (Prinn & Barshay 1977) and brown dwarfs such 
as GL 229b (Griffith & Yelle 1999). Diffusive tropo- 
sphcric mixing, in combination with detailed CO chem- 
istry, has recently been used to infer the water inventory 
in the deep Jovian atmosphere (Visscher ct al. 2010). 
Cooper & Showman (2006) parameterized the quench 
chemistry of CH4 in order to study its horizontal and 
vertical transport in their GCM of HD 189733b. The 
recent paper by Moses et al. (2011) discusses in de- 
tail the quench chemistry of H,C,N,0 molecular species 
in the relatively hot atmospheres of HD 189733b and 
HD 209458b. 

In our kinetics models we set thermochemical abun- 
dances as boundary conditions; these equilibrium abun- 
dance boundary conditions also define the metallicity 
of the system. We affix the 10 3 bar mixing ratios of 
the large carbon, oxygen and nitrogen reservoirs, CH4, 
H2O, CO, N2, and NH3, at their thermochemically de- 
rived values (here we are excluding sulfur), and set all 
other species to obey a zero flux condition at the lower 
boundary. The exact location of this lower boundary is 
unimportant, provided it is at depths much greater than 
the quench level (> 100 bars), and conditions (the high 
densities and temperatures) favor thermochemical equi- 
librium concentrations for practically all species. The 
nominal case has a solar abundance atmosphere (( = 1), 
vertical mixing with strength K zz = 1 x 10 8 cm 2 s _1 , 
and no photochemistry. In Figure 3 we compare an at- 
mosphere with vertical mixing to one purely in equilib- 
rium. Below 10s of bars, the mixing ratios converge, 
satisfying the condition that equilibrium concentrations 
have been reached kinetically. Now consider the abun- 
dances of quenched CO. At pressure levels deeper than 
10s of bars, the eddy mixing time, Ttrans, must be longer 
than the chemical loss timescale. As a check for internal 
consistency, we estimate 



Ttr 



L 2 

K zz 



8 x 10 5 s 



(3) 



where L is a fraction / of the scale height H, L = fH 
(Smith et al. 1998). We estimate / = 0.3 for both 
quenched CO and N 2 . To estimate T c h em for CO, we 
need to identify the rate-limiting reaction in CO and CH4 
interconversion. 



H + CO + M -> HCO + M #605 
H 2 +HCO^H 2 CO + H #234 



H + H 2 CO + M-> CH3O + M #611 

H 2 + CH3O -> CH3OH + H #351 

H + CH3OH -> CH 3 + H 2 #295 

H 2 + CH 3 -> CH 4 + H #61 

Net:3H 2 + CO^CH 4 + H 2 / 

This set of reactions is identical to the ones identified for 
CO quenching in Jupiter (Yung ct al. 1988; Visscher et 
al. 2010). The rate-limiting reaction is R351, the inverse 
of a hydrogen abstraction from methanol. The chemical 
loss timescale for CO is, 



T~chem,CO 



[CO] 



fc 3 5i[H 2 ][CH 3 0] 



(4) 



and h 



351 



where [X] denotes the concentration X, 

2.10 X l()-25 r 4.0 e -2470/T cm 3 mo i-l g -l ( Jo dkowski et 

al. 1999) the rate coefficient for R351. Figure 4 shows 
that equality of these two timescales, T c h em ,co ~ Ttrans, 
gives a CO quench-level of ~30 bars, which furthermore 
agrees well with quench-level depicted by the CO mixing 
ratio profiles in Figure 3. 

In an analogous manner the N 2 quench-level may be 
calculated by identifying the rate-limiting step in the se- 
ries of reactions that convert nitrogen to ammonia, and 
vice versa. These reactions are 



H + N 2 + M -> N 2 H + M #629 

H 2 + N 2 H -> N 2 H 2 + H #478 

H 2 + N 2 H 2 -> NH 2 + NH 2 #450 

2(H 2 + NH 2 -> NH 3 + H) #453 

Net : 4H 2 + N 2 -> 2NH 3 + 2H II 

In this N 2 ^ NH 3 sequence R450 is the rate-limiting 
step, involving the N abstraction from diazene, giving 
a timescale 



fN 2 



^~chem,N 2 



fc 45 o[H 2 ][N 2 H 2 ] 



(5) 



with reaction rate fc 450 = 2.06 x icr 07 T- - 93 e- 20614 / T , 
obtained from that of its reverse reaction (Stothard et al. 
1995). Calculating T c hem,N 2 above gives a N 2 quench- 
level of ^300 bars (see Figure 4), in agreement with 
the vertical profiles in Figure 3. The abovementioncd 
quench-levels for CO and N 2 are for the adopted eddy 
diffusion coefficient, K zz = 10 8 cm 2 s _1 . Increasing K zz 
to a very large value, 10 11 cm 2 s _1 , shortens the trans- 
port times considerably and increases the quench pres- 
sures of CO and N 2 to ~ 150 bars and <~ 620 bars, re- 
spectively. The effects of varying the quench level may 
be seen in Figure 2 - the atmospheric concentrations of 
the reservoir gases, CH4, H 2 and NH3, and quenched 
N 2 , are relatively insensitive to the location of quench 
pressure. However, varying the quench-level affects the 
concentration of CO and C0 2 by orders-of-magnitude. 
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Fig. 3. — Effects of vertical mixing on the distributions of H2O, CH4, NH3, N2, CO, CO2, and H. The dashed curves are the thermochemical 
equilibrium profiles for the f = 1 case from Figure 2 (middle panel). The solid curves are the vertical profiles derived kinetically with only 
eddy (K zz = 1 X 10 s cm 2 s _1 ) and molecular diffusion (no photochemistry) using the 1000 bar f = 1 mixing ratios as the lower boundary 
condition. Note that kinetically derived profiles begin to match the thermochemical equilibrium profiles at levels below a few 10's of bars. 
The rapid fall off of the solid curves near 1 /ibar is due to the sedimentation of the heavier molecules because of molecular diffusion. 



Vertical dredging of gases leaves a reasonably altered 
composition in the 1-0.001 bar region, the range of pres- 
sure levels wherein the infrared photosphere is located 
(e.g., Knutson et al. 2009; Swain et al. 2009). For ex- 
ample CO is up to a factor 10 4 more abundant than it 
would otherwise be. The deep quenching of N bearing 
gases causes NH3 to be surprisingly abundant, dominat- 
ing over the thermochemically favored N2. In contrast, 
the largest C and O reservoirs and optically the most 
active gases, CH4 and H2O, are largely unaffected. 

3.3. Photochemical Effects 

Photochemistry can significantly alter atmospheric 
composition in the upper portions. The combination of 
the ultraviolet flux and molecular absorption cross sec- 
tions gives the photolysis rates for all the species consid- 
ered here. The altitude of peak production/loss (in units 
of cm~ 3 s _1 ), set by the balance between the exponential 
fall-off of atmospheric density and the inward stellar UV 
attenuation, occurs near 1 /ibar (this is the well known 
Chapman function, see Yung & DeMore 1999 pg. 45). 
Primarily, photolysis breaks apart stable molecules into 
radicals, which can then react to alter the composition 
of the upper atmosphere. See Figures 5, 6 and 7 for the 
photochcmically derived mixing ratios. Table 1 compares 
the column mixing ratios from our models to the obser- 
vations over the 7 bar to 0.1 bar range probed by the 
observations. Figure 8 illustrates how photochemistry 
alters the upper atmosphere. The resultant mixing ratio 
profiles are compared with those obtained via thermo- 
chemical equilibrium (Figure 2), and by vertical mixing 



D 

rj2 



0. 



1 .0 



10.0 



I 00.0 



1000.0 











30 Bars 








300 Bars 







10 J 1 4 10 5 1 6 1 7 10 8 1 9 10 10 
Timescale [s] 

Fig. 4. — The blue curve is the CO chemical loss timescale 
calculated from Equation 4. The red curve is the N2 chemical loss 
timescale calculated from Equation 5. The dashed curve is the 
vertical mixing timescale from Equation 3 using a length scale of 
~ 0.3H estimated from the Smith et al. (1998) procedure.. The 
intersection of the vertical mixing timescale and the chemical loss 
timescale is the quench level for the given species as indicated by 
the horizontal dotted lines. 



(Figure 3). 



3.3.1. Atomic H & H 2 



Arguably, the most important radical in these atmo- 
spheres is atomic hydrogen. Its relatively large abun- 
dance (~75% above 1 /ibar, Figure 6) drives the bulk of 
disequilibrium chemistry in the upper atmosphere. As 
is seen in Figures 5-7, when the atomic H abundance 
increases with altitude, the concentration of disequilib- 
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Fig. 5. — Mixing ratios for important radicals (OH, NO, O, H, 
and CH3) that drive the photochemistry for three metallicities = 
O.l(top), ( = 1 (middle), and £ = 50 (bottom)). Note how the CH 3 
profile very nearly tracks the H profile because CH3 is a direct 
consequence of the oxidation of methane in R60 

rium species increases with it. Hydrogen attacks the 
large stable reservoirs, NH3 and CH4, to build these de- 
sequilibrium species. In the cold Solar System giants, 
atomic hydrogen is primarily produced by the photosen- 
sitized dissociation of H2 via heavier hydrocarbons, and 
the photodissociation of CH4 and ethylene C2H4. In hot- 
ter giant planets, as in GJ 436b, the atomic hydrogen is 
made primarily by the photodissociation of water (Liang 
et al, 2003, Line et al. 2010, Moses et al., 2010). This 
is because, unlike in the Solar System's giants, water 
is not sequestered in clouds and is readily available for 
photolysis. Its large UV cross section combined with a 
large thermochemical abundance, makes water the most 
important source of atomic hydrogen in GJ436b. The 
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Fig. 6. — Photochemically derived mixing ratios for the abundant 
species (H 2 0, CH 4 , NH 3 , N 2 , CO, C0 2 , and H) for £ = 0.1 (top), 
C = 1 (middle), and C = 50 (bottom). 

detailed mechanism for producing H is the photosensiti- 
zation of H 2 using water via, 

H 2 + \w -> OH + H R25 
H 2 + OH ->■ H 2 + H i?169 



Net : H 2 -> 2H III. 

This photosensitization is efficient because H 2 dissoci- 
ates out to —2000 A, whereas H 2 dissociates only out to 
-800 A. H 2 acts as a photon sink, with factor — 10 
more photons available for its photolysis, than for direct 
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Fig. 7. — Photochemically derived mixing ratios for the dise- 
quilibrium species, the hydrocarbons and hydrogen cyanide, for 
C = 0.1 (top), ( = 1 (middle), and £ = 50 (bottom) 

H 2 photolysis. Because ol these factors the net photo- 
sensitized destruction of H 2 by H 2 proceeds 5 orders- 
of-magnitude faster than the direct photolysis of H 2 , and 
3 orders-of-magnitude faster than the photosensitized de- 
struction of H 2 via the hydrocarbons. The mixing ratio 
of water itself is largely unaltered below 1/xbar levels. 

3.3.2. CH4, & Hydrocarbons 

Thermochemically, methane is the most abundant hy- 
drocarbon. Overall it is the fourth most abundant species 
after H 2 0, H 2 and He, and it is the parent molecule for 
the synthesis of all other hydrocarbons. Methane mix- 
ing ratios are > 10~ 4 at altitudes below the 0.1 mbar 
level, even for the lowest metallicities. The models gen- 
erally have methane mixing ratios at least 3 orders-of- 
magnitude higher than concentrations retrieved from the 



observations (Madhusudhan & Seager 2011). Although 
photolysis seems not to significantly modify methane 
abundances, it does produce large concentrations of the 
methyl radical, CH3; this radical is important in the syn- 
thesis of heavier hydrocarbons. CH 3 is formed by pho- 
tosensitized dissociation of methane. The free atomic 
hydrogen from scheme 177 readily attacks methane to 
produce H 2 and CH3. The trigger and pathway for this 
is: 



H 2 + hv -> OH + H R25 
H 2 + OH^H 2 + H i?169 
H + CH 4 -» CH 3 + H 2 R60 



Net : CH 4 ->■ CH 3 + H. IV 



The methyl radical's mixing ratios can be as high as 
~ 10~ 4 , as in the ( = 1 case (Figure 5). Due to the 
warmer upper atmosphere, relative to that in the so- 
lar system giants, the oxidation of methane (via R60) is 
more two orders-of-magnitude more efficient than direct 
photolysis. Because the forward reaction (R60) proceeds 
more sharply with rising temperature than the reverse 
(R61), hotter upper atmospheres (as in HD 189733b and 
HD 209448b) will have a tendency to destroy methane 
more readily, especially when there are large quantities 
of photochemically produced atomic hydrogen present. 
This photosensitized destruction of methane causes it 
to decline sharply above ~10 /Aars; this is well below 
the planetary homopause, but well above the infrared 
photosphere (Figure 8). It also drives the production of 
heavier hydrocarbons. Little to no heavier hydrocarbon 
(C„H m , where n,m > 2) is expected via vertical mixing 
alone, with mixing ratios remaining below ~ 10 -10 at al- 
titudes above 1 bar. Methane photosensitization (scheme 
IV) converts the carbon into ethylene (C 2 H4), acetylene 
(C 2 H 2 ), and ethane (C 2 H 6 ) via 



H 2 + hv 


->• OH + H 


i?25 


H 2 + OH 


-> H 2 + H 


i?169 


2(H + CH 4 


-»■ CH 3 + H 2 ) 


i?60 


CH 3 + CH 3 + 


M -> C 2 H 6 + M 


i?613 


H + C 2 H 6 - 


C 2 H 5 + H 2 


R70 


H + C 2 H 5 - 


-> C 2 H 4 + H 2 


R68 


H + C 2 H 4 - 


-> C 2 H 3 + H 2 


R85 


H + C 2 H 3 - 


-> C 2 H 2 + H 2 


i?64 



Net : 2CH 4 + 4H ->■ C 2 H 2 + 5H 2 V 

The net reaction ultimately produces C 2 H 2 , making it 
the most abundant heavy hydrocarbon. This scheme is 
different than the solar system gas giants where the most 
dominant pathway for producing acetylene involves the 
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binary collision between two 3 CH2 radicals. This differ- 
ence can again, be owed to the overwhelming abundance 
of atomic H from water photolysis which can readily re- 
duce the ethane produced R613 to acetylene. Over the 
range of metallicities considered (( = 0.1 to 50), the 
peak values of C 2 hydrocarbons occur between 10 and 
1 /ibars. These mixing ratios of C2H4, C2H2, C 2 H@ lie 
between 3 xl0" 7 -6xl0" 6 , 5xl0" 6 -4xl0" 4 , and5xl0~ 9 - 
6xl0~ 5 (Figure 7; for integrated columns see Table 1). 
For comparison, the peak values in Jupiter are, respec- 
tively, - 2 x 10~ 6 , 5 x 10~ 6 , and 20 x 10~ 6 (Moses et 
al. 2005). In the Solar System's giant planets, ethylene, 
acetylene, and ethane have strong mid-infrared strato- 
spheric emission features at 10.5, 13.7 and 12.1 ^tm re- 
spectively. These C2 species can lead to further synthesis 
of higher order hydrocarbons that can form hydrocar- 
bon aerosols (Zahnle et al. 2009). However, the vapor 
pressures for these species are high (many bars) at these 
temperatures, so it may be difficult to form such aerosols. 
Additionally, Moses et al. 1992 showed that supersatu- 
ration ratios of 10 to 1000s may be required in order to 
trigger condensation due to the lack of nucleation partic- 
ulates in Jovian type atmospheres, detectability 

3.3.3. CO & C0 2 

As described in §3.2, the CO abundance above 10 bars 
is determined by the reaction rate of scheme I, and the 
strength of vertical mixing. In the absence of incident 
stellar UV, a profile with a constant vertical mixing ratio 
up to the homopause is obtained. With incident UV 
radiation, there is a photochemical enhancement of CO 
near the 1 /zbar level, of up to a factor of 10 2 for the 
C = 1 case (Figure 6, 8). This high altitude enhancement 
is a property of the cooler atmosphere of GJ 436b; in hot 
Jupiter atmospheres, as in HD 189733b and HD 209458b, 
such enhancements or deficits will tend to be driven back 
towards equilibrium values. The carbon in this extra 
CO is ultimately derived from the CH4 reservoir, via the 
following reaction scheme: 



H 2 + his -> OH + H 


#25 


H 2 + OH^ H 2 + H 


#169 


H + CH 4 -> CH 3 + H 2 


#60 


H 2 + his -> O + 2H 


#26 


O + CH 3 -> H 2 CO + H 2 


#98 


H + H 2 CO^HCO + H 2 


#233 


H + HCO -> CO + H 2 


#213 



Net : H 2 + CH 4 -> CO + 2H 2 + 2H VI 

Scheme VI is driven by the water photolysis driven dis- 
sociation of CH4 to CH3 via scheme IV. Atomic O is 
produced by photolytic fragmentation of water (R26); 
the net absorption cross section for this branch is ~ 0.1 
that of the main branch in R25. The two radicals, O and 
CH3, form formaldehyde in R98, and followed thereafter 
by a two-step conversion to CO (R233 and R213). An 
enhancement of CO2 largely traces the enhancement of 
CO via: 



H 2 + his -> OH + H #25 
OH + CO^C0 2 +H #187 



Net : H 2 + CO -> C0 2 + 2H VII 

Photochcmically enhanced CO2 mixing ratios reach ~ 
10~ 4 at 1 /xbar for ( = 1. Column averaged mixing 
ratios arc 5 x 10~ 6 and 6 x 10~ 9 (see Table 1). This is 
low compared to the observed mixing ratios of ^lxl0~ 4 
and ^lxl0~ 7 , respectively. Increasing the metallicity 
to C = 50, increases the mixing ratios to ^lxl0~ 2 and 
~5xl0 -4 , suggesting that the observed CO and CO2 
columns are consistent with a metallicity enhanced to 
levels observed in Solar System's ice giant planets (Table 
!)• 

3.3.4. Nitrogen & HCN 

Ammonia and molecular nitrogen, N 2 , are thermo- 
chemically the two most stable species in a reducing 
atmosphere and their relative abundance within the 
1 — 0.001 bar pressure levels is dictated by quench chem- 
istry. Because it is relatively abundant, the addition of 
hot (quenched or otherwise) NH 3 (Tennyson et al. 2010) 
to the list of absorbers used for model fitting and retrieval 
may well be quite important. Other important N species 
are mainly photochemical byproducts, with HCN being 
the most abundant photochemically produced molecule 
between 1 and 0.1 mbar levels, having mixing ratios of 
typically 10~ 6 (C = 1) to 10~ 5 (C = 50) at 0.1 mbar. 
Peak HCN occurs well above the photosphcric levels, 
approaching 10~ 4 at 1 /ibar. The synthesis of HCN is 
initiated via water and ammonia photolysis, and com- 
pleted by subsequent reactions between the ammonia and 
methane derived radicals: 



H 2 + his OH + H #25 

H 2 + OH^H 2 + H #169 

H + CH 4 -> CH 3 + H 2 #60 

NH 3 + his -> NH 2 + H #43 

H + NH 2 -> NH + H 2 #455 

NH + CH 3 -I- CH 2 NH + H #685 

H + CH 2 NH -> H 2 CN + H 2 #655 

H + H 2 CN -> HCN + H 2 #663 



Net : CH 4 + NH 3 -> HCN + 3H 2 VIII 



We note that R43, the photolysis of ammonia to amino 
radical, is the most important pathway for NH 2 forma- 
tion at pressures greater than 10 ^ibar. At lower pressures 
this reaction is driven by ammonia photosensitzation, 

NH 3 + H NH 2 + H 2 , #454 

where the is H derived from H 2 photolysis. In con- 
clusion when water, ammonia and methane are present, 
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Fig. 8. — The effects of vertical mixing and photochemistry compared with the thermochemcial equilibrium profiles for methane, water, 
carbon monoxide, and carbon dioxide under solar abundances (C=l)- The dashed curves are the thermochemcially derived mixing ratios 
(Figure 2 middle panel). The solid curves are the mixing ratios with eddy mixing (as in Figure 3 middle panel). The dot-dashed curves 
are the photochemical mixing ratios (Figure 5 middle panel). Note that methane becomes photochemically depleted near 1 /ibar when 
compared to just vertical mixing (solid black curve). CO and CO2 are photochemically enhanced above the 10 mbar level. 



disequilibrium HCN is relatively abundant. The best 
chance for the detection of HCN is via the transmission 
spectroscopy of its vibrational fundamental bands at 3 
and 14 /im (Shabram et al. 2011). 

3.3.5. Sulfur 

Because atomic H attacks both CH4 and NH3, we ex- 
amine the role of H2S as a source of free H (Zahnle et 
al. 2009); S is isoelectronic with and similar in chemical 
properties to O, but has a considerably reduced primor- 
dial abundance, with S/O ~ 0.02. In a subset of models, 
we introduce the following (very restricted) set of sulfur 
reactions with accurate laboratory determined reaction 
rates: 



SH + H 2 -> H 2 S + H 



Net : H 2 -> 2H 



#712 



IX 



H 2 S 
H 2 S + H 
H 2 S + OH 



his 



■> SH 
SH + 
SH + 



f H 
H 2 
H 2 



#705 

#701, #702 
#703, #704 



H 2 S is an attractive source of free hydrogen due to 
its ability to photodissociate out to relative long wave- 
lengths, ^2600 A. It has a photolysis rate constant com- 
parable to that of H 2 0, and we find a 10 2 enhancement 
in H between the pressure levels of 1 bar and 0.1 mbar 
upon including these two sulfur species (Figure 9); the 
relevant reactions are: 



This enhanced H abundance is catalyzed by the photol- 
ysis of H 2 S (traced by the SH radical in figure 9, top 
panel). The atomic H reacts efficiently with CH4 in R60, 
producing an increased concentration of the radical CH 3 , 
which in turn drives hydrocarbon production (scheme V) 
near the 0.1 bar level. 

However, the free H in the middle atmosphere, does 
little to affect the CH4 mixing ratios; this is because the 
S/C abundance ratio is low. Sulfur would need to be en- 
riched by a substantial factor of ~20, over the solar S/C 
value, in order for H 2 S to have an appreciable impact on 
atmospheric CH4. Although the few considered sulfur 
species (H 2 S, SH) do not much impact the overall chem- 
istry, it is possible that another sulfur compound, such as 
SO, may act as a catalyst assisting in the conversion of 
reduced carbon into oxidized carbon. Previously, Moses 
(1996) has modeled the SL9 Jupiter impact and shown 
the importance of S in many reaction schemes involving 
both C and N species, and so the role of S chemistry in 
the hot extrasolar giants should continue to be investi- 
gated in the future (see Zahnle et al. 2009). 



H 2 S + \w -> SH + H 



#715 



4. DISCUSSION & CONCLUSIONS 
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TABLE 1 

Photochemical model abundances compared with the observations of Stevenson et al. 2010 (S10), Madhusudhan & Seager 
(2011) (MS10) and Beaulieu et al. (2010) (B10). The model abundances are the integrated column mixing ratios between 7 

BARS AND 0.1 BARS, THE PRESSURE LEVELS SAMPLED BY THE OBSERVATIONS, FOR O.lX, IX, AND 50X SOLAR ELEMENTAL ABUNDANCES. 



Molecule 


O.lx 




lx 




50 X 




MS10 


S10 


B10 


CH 4 


7.66 xlO" 


-05 


7.90x10" 


-04 


2.96 xlO 


-02 


(3-6)xl0-° 6 


lxl0~ 07 


5 xlO" 04 


CO 


4.22x10" 


-08 


4.29x10" 


-06 


8.56x10" 


-03 


(3-100) xlO" 08 


(l-7)xl0-° 4 




C0 2 


7.74x10" 


-12 


6.09x10" 


-09 


5.44x10- 


-04 


(1-10) xlO" 07 


(1-10) xlO- 07 




H 2 


1.25x10" 


-04 


1.26x10" 


-03 


5.09x10" 


-02 


<lxl0" 03 


(3-100) xlO" 06 




HCN 


4.84x10" 


-10 


3.09x10" 


-08 


8.41 xlO 


-06 








C2H2 


1.21x10" 


-14 


1.18x10" 


-12 


2.10 xlO 


-09 








NH 3 


1.45x10" 


-05 


1.06x10" 


-04 


6.54x10" 


-04 








H 2 S 






3.22x10" 


-05 













We have developed a ID "thermochemical and photo- 
chemical kinetics with transport" model following Viss- 
cher et al. (2010) and recently, Moses et al. (2011) for 
extrasolar planet atmospheres. We use a compilation 
of bidirectional reactions of the five most abundant ele- 
ments to model both the equilibrium and disequilibrium 
portions of the atmosphere. Using detailed balance with 
both forward and reverse reactions, allows our model 
to reach thermochemical equilibrium kinetically, thereby 
obviating the need to choose ad hoc lower boundaries for 
multiple quenched species, and allowing a seamless tran- 
sition between the transport dominated and the chemical 
equilibrium zones. A limitation is that we adopt a static 
temperature structure; a future improvement would al- 
low the iterative adjustment and co-evolution of the tem- 
perature structure with the chemistry. Also, the eddy 
diffusivity profile K zz (z) is poorly constrained, and is 
essentially a free parameter in any of these models. 

We have applied our models to study the atmosphere 
of the transiting Neptune-like planet GJ 436b. The ele- 
mental abundance of atmosphere, a key input parameter, 
is relatively uncertain, but mass-radius constraints sug- 
gest that GJ 436b must be enriched to at least 10 x solar 
levels. We model a range of atmospheric enrichment to 
cover this instrinsic uncertainty; we observe the trends 
when varying £, and rule out the possibility that inter- 
mediate values of £ would spring any surprises. The UV 
fluxes of stars other than the Sun are often difficult to 
obtain. M-dwarf hosts can be chromospherically hyper- 
active, and because UV photolysis may drive the deple- 
tion of weakly bonded molecules such as CH4, NH 3 and 
H2S, it is important to have an accurate UV estimate 
for GJ 436. We use a combination of GALEX and HST 
UV fluxes along with Rosat and XMM-Newton soft X- 
ray fluxes to bound the UV continuum and line emission 
of GJ 436. 

The GJ 436b model atmospheres show that a combina- 
tion of photochemistry, chemical kinetics and transport- 
induced quenching drives the composition well out of 
equilibrium. While equilibrium conditions are main- 
tained in the deep, hot, troposphere (below a 10s of bars 
for CO ^ CH 4 , and 100s of bars for N 2 ^ NH 3 ), the 
composition of the middle atmosphere is altered by the 
dredging up of quenched gases such as CO and NH 3 . 
The effects of transport disequilibrium are prominent in 
cooler planets such as GJ 436b because the quench points 
for major species depend on the temperature. As it 
gets colder, the pressure points for quenching are pushed 



deeper into the atmosphere due to the longer intercon- 
version timescales from one species reservoir to another. 
In contrast to the quenched species (CO, CO2, NH 3 ), 
the effect of vertical mixing on the reservoir gases such 
as CH4 and H 2 is relatively feeble. 

The reservoir gases H 2 and CH4, and NH 3 are largely 
unaffected by photochemistry because of their (a) large 
abundances, and (b) rapid recycling. Nevertheless, it is 
their photolysis that drives the bulk of the disequilibrium 
chemistry in the upper atmosphere producing CH4 and 
NH 3 sinks such as heavier hydrocarbons (such as C2H2, 
etc.) and simple nitriles (such as HCN). Much as in the 
hot Jupiters (Liang et al. 2003), H is the most impor- 
tant and active atom in the bulk of the atmosphere; it 
is created by the photosensitized destruction of H2 , cat- 
alyzed by the presence of H 2 and H 2 S. The latter gas, 
though less abundant than water, is important because of 
its ability to capture incident starlight photons out wave- 
lengths as long as 2600 A. In most models, H replaces H2 
as the most abundant species in the atmosphere above 
the planetary homopause at ~ 1 /ibar. Because CH4 is 
the largest C carrier in the planet's UV photosphere, we 
create abundant C 2 compounds (Figure 7) despite the 
relatively efficient hydrogenation back to CU4. Species 
such as acetylene, C 2 H 2 , formed in abundance in our en- 
riched models, are precursors for potential hydrocarbon 
soot formation in the upper atmosphere (as opposed to 
the hotter Jupiters such as HD 209458b and HD 189733b, 
wherein CO carries the bulk of carbon in the strato- 
sphere). Our reaction lists for hyrdocarbon chemistry 
are truncated at C 2 , and so we do not synthesize C 3 and 
heavier hydrocarbons and nitriles explicitly. 

Within the range physical and chemical processes cap- 
tured in our models, and the considered reaction sets and 
their kinetics, we find it difficult explain the observations 
suggesting a methane-poor GJ 436b. Except above 1 
/xbar pressure levels where CH4 is photochcmically con- 
verted to CO, HCN and C 2 hydrocarbons, it remains the 
predominant C reservoir in the lower atmosphere and in 
the region of the IR photosphere. The observed abun- 
dances of quenched CO and C0 2 are in agreement with 
an atmosphere enriched to levels intermediate between 1 
to 50 times solar (as in Madhusudhan & Seager 2011). 
The depleted water may either contrarily suggest a sub 
solar metallicity (Table 1), or skewed heavy metals ra- 
tios; the latter is a possibility which we have not con- 
sidered herein as there are far too many combinations to 
explore. In the 1 x solar models, the methane abundance 
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Fig. 9. — Photochemically derived mixing ratios for f = 1 in 
the presence of sulfur species (H2S and HS) for the radicals (top), 
abundant species (middle), and disequilibrium species (bottom). 
Compare this Figure to the ( = 1 cases in Figures 5, 6, and 7 to 
see the effects of H2S on the mixing ratios. Note that the abun- 
dance of atomic H is enhanced by orders of magnitude between 1 
and 10 — 4 bars as a result of scheme IX. This H increase enhances 
the hydrocarbon abundances significantly in this portion of the 
atmosphere. 

is consistent with the values retrieved by Beaulieu et al. 
(2010) (Table 1) using transit observations. We suppose 
it is possible that a more complete inclusion of other rel- 
atively abundant elements such as S and P, or distorted 
elemental ratios (C/O or O/S etc.), or ill-understood 
chemistry and exotic processes (not considered herein, 
such as the 3 dimensionality of the problem) could do 
more to explain the chemistry of this enigmatic atmo- 
sphere. 

We agree with Moses et al. (2011) that quench level 
arguments can be used to predict abundances, so long 
as this is done with the appropriate level of caution. By 
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Fig. 10. — Mixing ratios for, CO, CH4, CO2, and H2O assuming 
a T=1200 K isothermal atmosphere, f = 5, and K zz = l X 10 
cm 2 s _1 , with no photochemistry. The observed methane fall-off is 
due to high temperatures alone; high temperatures imply a short 
chemical loss time for of CH4. Because of the large transport time, 
CH4 and CO are nearly in equilibrium at all altitudes. 

this, we mean that the relevant rate-limiting reaction 
must necessarily be identified in order to properly cal- 
culate the timescale for chemical loss. Also, quenched 
gases do not share a common quench level and assuming 
so can result in gross under- or overestimation of their 
abundances. For example, as shown herein, N2 and CO 
have vastly different quench levels. For the moderate 
to high levels of incident UV flux, photolysis generates 
high concentrations of secondary byproducts, but does 
not significantly alter the abundances of the reservoir 
gases; in our estimation photochemistry cannot alter the 
dayside methane budget. Hotter atmospheres with slug- 
gish vertical mixing and hot stratospheres are required 
for severe methane depletion. For example, in figure 10, 
we approximate such as atmosphere as isothermal with 
T = 1200 K, C = 5, and K zz =l x 10 6 cm 2 s~\ and 
with zero UV irradiation (similar to models by Zahnle 
et al. 2009). In this hypothetical atmosphere there is 
relatively little quenched methane. At T = 1200 K and 
low pressures, the rate determining step for CH 4 — > CO 
(reverse of R351) is faster than the vertical transport 
time throughout the atmosphere, allowing the CH4 to 
be in thermochemical equilibrium with CO everywhere 
(Figure 10). Since equilibrium conditions apply, the P 2 
term in equation 2 results in the rapid vertical fall-off of 
CH 4 . 

The models presented herein are by no means re- 
stricted in applicability to GJ 436b like Neptunes, and 
much of the modeled chemical state may be generalized 
to H/He dominated planets in the 500-1000 K temper- 
ature range. In this regime CH4 is the primary carbon 
carrier and CO is quenched. The reverse is true in hotter 
atmospheres, T > 10 3 K, where CO is the primary car- 
bon carrier and CH4 is quenched. NH3 is quenched deep 
in the atmosphere and can be quite abundant in the pho- 
tosphere. Higher hydrocarbons and HCN are produced 
photochemically in relatively high abundances at mbar 
to /ibar pressures. Similarly, an enhancement of CO and 
CO2 over the quench concentrations, driven by the pho- 
tolysis of H2O, is observed in the high atmosphere. Wa- 
ter is in gaseous phase and abundant, and not condensed 
out as it would be in cooler atmospheres. GJ 1214b, a 
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T ~ 500 K low super Earth or mini Neptune, also orbit- 
ing an M dwarf primary (Charbonneau et al. 2009; Sada 
et al. 2010; Bean et al. 2010; Carter et al. 2011; Desert 
et al. 2011), falls in this regime of warm atmospheres. 
If GJ 1214b is in possession of a reducing H-He atmo- 
sphere (Croll et al. 2011; Crosshcld et al. 2011), much 
of the atmospheric chemistry would be analogous to that 
in GJ 426b; this, however, is speculative as there is much 
current debate over the bulk composition of GJ 1214b. 
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